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Abstract 

We prove the existence and stability of flat steady states of the 
Vlasov-Poisson system, which in astrophysics are used as models of 
disk-like galaxies. We follow the variational approach developed by 
Guo and Rein fS", "B", T for this type of problems and extend previous 
results of Rein JA^. In particular, we employ a reduction procedure 
which relates the stability problem for the Vlasov-Poisson system to 
the analogous question for the Euler-Poisson system. 

1 Introduction 

In astrophysics, galaxies or globular clusters are often modeled as a large 
ensemble of particles (stars) interacting only by the gravitational field which 
they create collectively. In such systems collisions among particles are sufh- 
ciently rare to be neglected. Hence in a nonrelativistic setting the particles 
move on trajectories determined by Newton's equations of motion 

X = V, V = -VxU{t,X), 

where U{t,X) denotes the gravitational potential of the ensemble, tGM is 
time, and X,VgM'^ denote position and velocity; it is assumed that all 
particles have the same mass which is normalized to unity. To describe the 
time evolution of the ensemble, the density function F:RxM^xR^^R^ 
on phase space is used. It obeys a conservation law known as the Vlasov 
equation, 

^ + V-VxF-VxU-VvF = 0, 
ot 
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which can be understood as an incompressibihty condition of the " fluid" in 
phase space, also known as Liouville's theorem. The gravitational potential 
U is induced by the spatial density 



R{t,X):= [ F{t,X,V)dV 



via Newton's law of gravity 




f R{t,Y) 
V^\X-Y\ 



dY. 



In the present paper we are interested in a situation where extremely flat- 
tened objects such flat galaxies are to be modeled. We therefore assume 
that all particles are concentrated in a plane, say the (xi,a;2)-plane, with 
velocity vectors tangent to it. If this holds initially and if the only force 
acting on the particles is their mutual gravitational attraction, then the 
particles stay in that plane, and we can introduce a new, flat particle den- 
sity / :]R X X R2 which is related to the density on the full phase 
space through 



Here 6 denotes the Dirac distribution, and X = (xjX^), V= (vjV^) with x,vG 
M^. However, the particles still interact by the three dimensional Newtonian 
gravitational potential. Therefore, the Vlasov-Poisson system for / reads 



where x,uGM^. We refer to this system as the flat Vlasov-Poisson sys- 
tem. From the mathematics point of view it consists of the two dimensional 
Vlasov equation coupled with the l/|3;|-type potential typical for three space 
dimensions. Since the l/|x| -singularity is integrated only over this kind of 
coupling makes the system more singular and mathematically more difficult 
to analyze than the regular, three dimensional system. 

The aim of the present investigation is to prove the existence of a large 
class of non-linearly stable steady states of the flat Vlasov-Poisson system. 
To do so we follow the approach developed by Guo and Rein El EI 



F{t,X,V) = f{t,x,v)6(x3)6{v3). 



-£ + v-V,f-V,U-VJ = 0. 
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in the regular, three dimensional situation. We prove that under suitable 
assumptions on a prescribed function <I> : [0,oo[^ [0,oo[ the energy-Casimir 
functional 

Wc(/) ^^jj \v\'^f{x,v)dvdx-^ jj jj '^^^1^'*/^^'"'^ dvdxdwdy 

+ jj ^{f{x,v))dvdx 

has a minimizer /o subject to the constraint 

jj f{x,v)dvdx = M, 

where M > 0, the total mass of the resulting steady state, is prescribed. The 
Euler-Lagrange relation for this variational problem implies that 

Ux,v) = <i,{E). 

Here the particle energy E is defined as 

E{x,v) = ]^\v\^ + Uq{x) 

with Uq the potential induced by /o, and the function (j) is determined by <I> 
and a Lagrange multiplier. The point now is that for the time-independent 
potential Uq the particle energy E and hence also (l){E) is constant along 
particle trajectories and hence a solution of the time-independent Vlasov 
equation. Hence /o is a steady state of the flat Vlasov-Poisson system. The 
fact that /o minimizes the energy-Casimir functional Tic can then be used 
to derive a non-linear stability property for this steady state. 

In this approach has already been used to construct stable steady 
states of the flat Vlasov-Poisson system. In the present paper we obtain 
a number of improvements and extensions of this earlier result. Firstly, 
we use a reduction procedure for proving the existence of a minimizer of 
Tic- This approach is mathematically more elegant and adequate, since 
the reduced functional lives on the set of spatial densities p, and the main 
difficulty in the variational problem lies in the potential energy part which 
does not really depend on / but only on the spatial density induced by /. 
More importantly, the reduced variational problem is of interest in its own 
right since it provides a stability result for the flat Euler-Poisson system 
which is the fluid dynamical analogue of the kinetic Vlasov-Poisson system. 
For the reduction procedure to work the function $ has to satisfy certain 
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growth conditions. An example of a steady state which violates this growth 
condition is the so-called Kuzmin disk which is known in the astrophysics 
literature and was not covered by previous results. The Kuzmin disk will be 
investigated in a companion paper Secondly, in the perturbations 
admissible in the stability result had to be supported on the plane and in 
addition had to be spherically symmetric. In the present paper we remove 
the latter, unphysical restriction. It is desirable to remove also the restriction 
that the perturbations have to live in the plane, but that is much harder 
and is still under investigation. Lastly, we relax the assumptions on $ the 
main one being that $ be strictly convex so that we cover a larger class of 
steady states, and we obtain stability estimates in stronger norms than were 
obtained previously. 

The paper proceeds as follows. In the next section we introduce vari- 
ous functionals and the reduced version of the variational problem, and we 
establish the connection between the original and the reduced variational 
problem. In the third section we establish the existence of a minimizer 
to the reduced problem using a concentration-compactness argument; no- 
tice that the variational problem — both reduced and original — is non-trivial 
since the energy-Casimir functional is not convex and is defined on func- 
tions supported on M'^ or respectively. In Section 4 we derive our sta- 
bility result, and in the final section we consider the stability result for the 
Euler-Poisson system which arises from the reduced functional. 

2 Energy-Casimir functionals and reduction 

For p = p{x) measurable we define the induced gravitational potential and 
potential energy as 

J \x-y\ 

in the sequel integrals J without a subscript always extend over M?. It 
will also be useful to introduce the bilinear form which corresponds to the 
potential energy, i.e., for p,cj:M^^]R measurable, 

1 ffp{x)a{y) 
Mpot:=^jj ^^dydx, 

so that in particular Epot{p) = —{p,p)pot- For the convenience of the reader 
we collect the main estimates for potentials, potential energies, and the 



4 



above bilinear form, which we will need. 

Lemma 2.1 If L^^^{M?), then UpG L^{M?), and there exists a constant 
C>0 such that for all yOGL^/^(R^) the estimates 

||C/p||4<C||p||4/3, -E^M<C\\p\\l^, 

hold. The bilinear form (•,-)pot defines a scalar product on L^/^(M^) with 
induced norm 

\\p\\:>ot:={p,p)lil = {-E^ot{p))'/', 

in particular, 

(/5,o-)pot < {Epot{p)Epot{(y)y^'^ = ll/ollpotlkllpot- 
Proof. Since l/\-\eLl{R^), the weak L2 space, the assertions on Up follow 
by the generalized Young's inequality |H| 4.3]. The estimate for the potential 
energy is nothing but the Hardy-Littlewood-Sobolev inequality [HI 4.3] and 
follows by Holder's inequality, and so does the fact that (•,-)pot is defined on 
L'^^^{M'^). The positive definiteness of (•,-)pot can be shown exactly like the 
positivity of the Coulomb energy in the three dimensional case, cf. 9.8]. 

□ 

Let f = f{x,v) be a measurable function on phase space. We define the 
induced spatial density, gravitational potential, and potential energy as 

Pfix)-= j fix,v)dv, Uf:=Upj, Epot{f):=Epot{pf)- 

In addition, we define the kinetic energy 



i^kin(/) : = 




f{x,v) dvdx, 



the so-called Casimir functional 




(^{f{x,v))d.vdx 



with $ : [0,00 [0,oo[ prescribed, and the energy-Casimir functional 

nc{f):=C{f) + E^^{f) + EpM- 

The total energy E\^\^ + Epot as well as the Casimir functional C and hence 
also their sum Tic are conserved along sufficiently regular solutions of the 
flat Vlasov-Poisson system. As regards we assume for the moment that 

$gC^([0,oo[) is strictly convex, ^>(0) = $'(0) = 0, lim <l>(r/)/r/ = oo. 



5 



These assumptions make ^ non-negative and <!>' a bijection on [0,oo[. 

Our aim is to show that the energy-Casimir functional TCq has a mini- 
mize! in the constraint set 

^M:=|/G4(M^)|^kin(/)+C(/)<oo, pfeL^/\R^),JJ / = m|, 

where M>0 is prescribed, and the subscript + indicates that only non- 
negative functions are considered. Since the troublesome term in the func- 
tional is the potential energy which actually depends only on the spatial 
density induced by / we introduce a reduced variational problem for a func- 
tional which is defined in terms of spatial densities p. For r > we define 



g{v)dv-- 



and 

^'(r):=inf / ( -{v]"^ g{v) +^{g{v))] dv. 
g&QrJ \2 J 

The idea behind this construction is to first minimize the energy-Casimir 
functional over all functions f{x,v) which upon integration in v give the 
same spatial density p, and then minimize with respect to the latter in a 
second (and main) step. This approach was introduced in |121ll7j . 

The reduced variational problem is to minimize the reduced functional 



nUp):= J ^ip{x))dx + EpM 



over the set 



jr^:=|pGL^/3nL^(M2)| j {p{x)) dx < oo , J p{x)dx = M 

We need to establish a relation between minimizers of the original functional 
and minimizers of the reduced one. Here we can essentially follow the cor- 
responding results proven for the three dimensional case in [El- First of all 
we explore the relation between $ and ^. For a function /i : M^] — 0x3,00] 
we denote by 

r(A):=sup(Ar-/i(r)) 

reM 

its Legendre transform. In what follows constants denoted by C are always 
positive, may depend on ^ and M, and may change their value from line to 
line. 
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Lemma 2.2 Let $ and ^ be as specified respectively defined above, and 
extend both functions by +00 to the interval ] — oo,0]. Then the following 
holds: 



and in particular $* (A) = = ^'*(A) for all A < 0. 

(b) ^'eC^([0,oo[) is strictly convex, and ^(0) = ^''(0) = 0. 

(c) Let k>0 and n = k + l. 

(i) If^{f) = Cf^+^/''forf>(}, then ^{p) = Cp^+'^/'' for p>0. 

(ii) // $(/) > C/i+V'^ forf>0 large, then ^{p) > Cp^'^^l'^ for p>0 



(iii) // $(/) < C/i+Vfc forf>0 small, then ^{p) < Cp^+^l'' for p>0 



If the restriction to large or small values of f can be dropped then the 
corresponding restriction for p can be dropped as well. 

Proof. By definition 



(a) ForXeR, 




large. 



small. 




1 



2 



) 



dv. 



For A< both sides of this estimate arc zero, so consider A > 0. If l^l > \/2A 
then supy>Q [•••]= and for < \/2A the supremum is attained at y = yv '■= 
($')~^(A— P). Hence with the definition 
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we obtain the reversed estimate, and part (a) is established. Part (b) is 
standard for Legendre transforms, and we refer to jl2l Lemma 2.2] for the 
details. As to (c), if we assume that <l>{f) > Cf^+^'^ for / > large, we find 
that $(/) > Cf^+^/^ - C for / > 0. Hence for A > 0, 

(A) = sup (A/ - $(/)) < C + sup (a/ - C7/i+i ) < C + C A'=+^ , 

/>0 />0 ^ ^ 

and 



^*{X)= [ <^>*(x-hvAdv<C [ i+(x-l\vA 

J\v\<V2X \ ^ J J\v\<V2\ \ ^ J 



<CX + C I (A--|t>|2) dv<C + CX''+^ = C + CX^+''. 



\v\<V2X 

Using the fact that ^** = ^' we obtain the estimate 

^ (p) = sup(pA -^*{X))>-C' + sup(/)A - (:7A^+") = Cp^+^Z" - C, 

A>0 A>0 

which proves (c) (ii). The remaining estimates are shown in a similar way. 

□ 

The relation between the minimizers of Tic and Ti^ is as follows. 
Theorem 2.3 (a) For every function f £ J^m , 

nc{f)>nh{pf), 

with equality if f is a minimizer of Tic over Tm ■ 

(b) Let poG-^/ « minimizer of TCq with induced potential Uq. Then 
there exists a Lagrange multiplier EqGM. such that the identity 



Po- 



i^')-\Eo-Uo) ,Uo<Eo 
,Uo>Eo 



holds almost everywhere. The function 

/o:=i ^ ^ \ ''^J'wtthE = E{x,v):=-\v\^ + Uo{x) 



,E>Eo ^ ' ^ 2 

is a minimizer of Tic 
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Proof. Since the proof follows the same lines as jl21 Thm 2.1] we only 
indicate the main arguments. The estimate in part (a) follows directly from 
the definitions. Next one can show that if f gJ^m is such that up to sets of 
measure zero, 

$'(/) = ^o-^>0 where />0, and - < where / = 0. (2.1) 

with E defined as in (b) but with Uf instead of Uq and Eq a constant, then 
equality holds in part (a). If / is a minimizer of Tic, then the Euler-Lagrange 
equation implies that / is of the above form for some Lagrange multiplier 
£"0, and equality holds in (a). The relation of po and Uq in part (b) is nothing 
but the Euler-Lagrange equation for the reduced variational problem. If /o 
is defined as in (b) then po = P/o) particular, foGJ^M, and 1)2. 1() holds by 
definition of /q. Hence equality holds in (a) for /o so that by part (a) for 
any other / G J^m, 

nc{f)>nh{pf)>nHpo)=nc{fo), 

which means that /q minimizes TCq- □ 

Remark, (a) In the next section we show that under suitable assumptions 
on ^ which can be translated into corresponding assumptions on <I> the 
reduced variational problem has a solution pQ. The minimizer /o obtained by 
the lifting procedure in part (b) of the theorem depends only on the particle 
energy E. The latter is for the time-independent potential Uq constant along 
characteristics of the Vlasov equation, and hence /o is a steady state of the 
flat Vlasov-Poisson system. 

(b) If TCq has at least one minimizer in and if /o G J-m is a minimizer of 
Tic, then one can show that pq:= pf^G J-\^ is a minimizer of Tf^. This map 
is one-to-one between the sets of minimizers of Tic in and Ti^Q in 
and is inverse to the mapping pQ 1— > /o described in part (b) of the theorem. 

3 Existence of a solution to the reduced varia- 
tional problem 

In the present section we prove that the reduced energy- Casimir functional 
Ti^Q has a minimizer in the constraint set 



^{p{x))dx <Qo, / p{x)dx = M 



where M>0 is prescribed and "if satisfies the assumptions ^' G C^([0,oo[), 
5'(0) = ^''(0)=0 and 
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(^'1) ^' is strictly convex, 

(^2) ^(p) > Cpi+V" for p> large, 

(^'3) ^'(p) < Cpi+i/"' for p>0 small, 

with growth rates n,n'E]0,2[. The core of the proof is a concentration- 
compactness argument to show that along a minimizing sequence the matter 
cannot spread out but has to remain concentrated in a finite region of space. 
First however we show that the energy-Casimir functional is bounded from 
below in such a way that minimizing sequences are bounded in a suitable 
LP space. 

Lemma 3.1 Under the above assumptions on ^ and for p£j^, 
j p^+^/"dx<C7 + C j ^{p)dx, 

and 

n/2 



mp)> j ^{p)dx-c-c(^j ^{p)dx 



In particular, 

^M:=infWc>-oo. 

•'M 

Proof. The first estimate follows by assumption (^'2) and the fact that 
J p = M. By Lemma l2. II and interpolation, 

J 'f{p)dxj , 

and since < n < 2 the proof is complete. □ 
We note an immediate corollary. 

Corollary 3.2 Any minimizing sequence of TCq in J-J^ is bounded in 
2^i+i/n^]^2^ and therefore contains a subsequence which converges weakly 
in that space. 

The concentration-compactness argument mentioned above relies on the be- 
havior of 'Hq{p) if p is scaled or split into several parts. We start with the 
latter; in the sequel Br denotes the open ball of radius R>0 about the 
origin. 
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Lemma 3.3 Let pEJ-]^^. Then for R>1, 

sup / p{x) dx>^ (-2Spot (p) - M'R-' -C\\p\ ) . 
Proof. We split the potential energy as follows: 

JJ\x-y\<l/R \X-y\ JJl/R<\x-y\<R JJ\x-y\>R 

=: /1+/2 + /3. 

By Holder's and Young's inequalities we obtain estimates 

/l< ||p||l+l/n||P*(lBl/Kl/|-|)l|n+l < l|p||?+l/„l|lSi/«V|-|ll(n+l)/2 

<c\\p\\i^,/^R-^'--y^-+'\ 

h^R p{x)p{y)dxdy = MR sup / p{x)dx, 

JJ\x-y\<R a&?Ja+Bii 
I3<M^R~^. 

We insert these estimates into the formula for — 2£^pot and rearrange terms 
to obtain the assertion. □ 

Next we investigate the behavior of the reduced functional under scalings. 
Lemma 3.4 (a) For cAiery M >0, /i^<0. 

(a) For every 0<Id<M the estimate h^> {MjMfl'^W^ holds. 
Proof. For peJ^ and a,6>0 we define p{x) :=ap{bx). Then 

pdx = ab~'^ / pdx, 



J *(p)dx = b""^ J *(a/9)dx. 

To prove part (a) we fix a bounded and compactly supported function p G 
and choose a = b^ so that pGj^j^ as well. By (^3) and since 2/n'> 1, 

Wc(p) = 6~'/ *(6V)dx + 6Epot(p) < Cfc^/-' + 6Epot(p) < 
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for b sufficiently small, and part (a) is established. As to part (b), we take 
a = l and 6= (M/M)i/^ > 1. For this choice of parameters the mapping 
3 p ^ p £ J-jY is one-to-one and onto, and the estimate 

>b-^y ^(p)dx+E^M) = (^^) ' ni{p) 

proves the assertion of part (b). □ 



Corollary 3.5 Let {pi)cJ-]^,j be a minimizing sequence o/TCq. Then there 
exist 5o > 0, Ro> 0, and a sequence of shift vectors (aj) C such that for i 
sufficiently large, 

/ Pi{x)dx>6o. 
Proof. By Corollary 13.21 (| |i+i/n) is bounded. By Lemma 13.41 fa). 

EpoM<ni{p^)<h\,<o 

for i sufficiently large, and the assertion follows by Lemma 13.31 □ 

This corollary only shows that along a minimizing sequence not all matter 
can spread uniformly. In the proof of the existence theorem below we shall 
actually see that the matter remains within a ball of finite radius up to 
spatial shifts and an arbitrarily small remainder. In such a situation we 
have the following compactness result: 

Lemma 3.6 Let [pi)<zT\.j he such that 

p^^Po weakly in l1+i/"(M2) 
and such that the following concentration property holds: 

Ve>03i?>0: limsup / pi{x)dx < e. 

i^oo J\x\>R 

Then 

EpotiPi- po)^0 and Epot{pi) ^ Epot{po), i^oo. 
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Proof. By weak convergence po>0 a. e., and / po<M. We define o"i:= 
Pi — po so that cTj ^ weakly in L^+^/" 0^"^), the concentration property holds 
for |(Tj| as well, and / |(Tj| < 2M. We need to prove that 

JJ \x-y\ 
which is the first assertion. Since 

Epot{Pi)-Epotipo) = Epotipi- Pq) + J Ufy^{pi-po), 

the fact that ?7p(, G L^(M?) together with the weak convergence of pi implies 
the second assertion. For S>0 and i? > we split the domain of integration 
into three subsets defined by 

\x-y\ <5, 

\x-y\ >SA{\x\>R\/\y\>R), 
\x-y\>d A\x\ < RA\y\ < R, 

and wc denote the corresponding contributions to /j by /j^i, Ii^2, Ii,3- Young's 
inequality implies that 



l-^i.ll < C\ \cTi\\l^yJ\lBs I • I ^ I |(n+l)/2 < C 



2/(n+l) 



which can be made as small as we wish, uniformly in i and R>0, by making 
(5 > small. For d>0 now fixed. 



Mi,2| <— ^ / \ai{x)\dx 

^ J\x\>R 



\x\>R 

which becomes small for i^oo by the concentration assumption, if we 
choose R>0 accordingly. Finally by Holder's inequality, 



-^1,3 1 = 



J ai{x)hi{x) 



dx 



< 



kj||l+l/nll^i||l+n <C'||^i||l+n, 



where in a pointwise sense, 

J\x-y\>5 \X-y\ 

Awe to the weak convergence of cjj and the fact that the test function against 
which (Tj is integrated here is in L^^". Since \hi\< ^1^^ uniformly in i. 
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Lebesgue's dominated convergence theorem implies that /ij — > in L^^", and 
the proof is complete. □ 

We have now assembled all the tools we need to prove the existence of a 
minimizer of the reduced functional. 

Theorem 3.7 Let [pi)(ZJ-\,.j he a minimizing sequence o/W^. Then there 
exists a sequence of shift vectors {a-i) CM^ and a subsequence, again denoted 
by (pi), such that for every e>0 there exist R>0 with 

/ pi{x)dx>M-e, ieN, 

Jai+BR 

Ta^pr.=Pi{- + a^)^po weakly in l1+i/"(r2), i^oo, 

/ po{x)dx>M-e. 
Jbr 

Finally, 

Epot{TaiPi- po)^0, 
and pq G J^j is a minimizer of TCq . 

Proof. We split p G as follows: 

P='^BR^P+'^BR^\BR^P+'^R-i\BR.^P = - Pl+ P2 + p-i- 

The parameters Ri < R2 of the split are yet to be determined. Recalling the 
definition of the bilinear form (•,-)pot) 

(P) = iPl )+^c{p2)+ nh (ps) -2{pi+ P3, P2)pot - 2(pi , p3>pot • 

If we choose R2 > 2Ri , then 

C 

(/5l,P3)pot<^. 

By Lemma I2. II and interpolation, 

{Pl + P3,P2)pot < ||/Ol+P3||pot||/02||pot 

< C'||/9l + /)3||4/3l|P2||pot <C||p||["'^/J/''||p2||pot- 



If we define 

Ml. 



J pi{x)dx, 1 = 1,2,3, 
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then Lemma 13.41 (h) and the estimates above imply that 



M J \M J \M J 

(n+l)/. 
1+1/n 

1+1/n 

< Chl,M^M^ + c 1 + 1 IpI I S";^)^/^ 1 1/,2 1 Ipot 
Here we used that for some constant C > the following inequality holds: 



+c(i?2^+iiHi;:trii^2iu 

c 

< ^(MiMs+MsMg + MiMg)/! 



3/2 



+ + < 1 — C(xy + x2; + yz) for x,y,z>0 with x + y + z = 1. 



Now we consider a minimizing sequence {pi)<zJ-\j of WJj and choose shift 
vectors (aj) C M^, (5o > 0, and i?o > according to Cor. 13.51 Since all our func- 
tionals are invariant under spatial translations the sequence T^-pj = Pi{-+ai) 
is again minimizing and hence bounded in L^+^/"(]R^) so that up to a sub- 
sequence we can assume that it converges weakly to some p^^ L^'^^^'^i^). 
We choose Ri > Rq so that by Cor. 13. 5( Mj^i > 5o for i large, and 



-Ch'MSoMi^S < Ci?2 ^ + C| |po,2| Ipot + C\ \pi,2 - P0,2\\pot + rL'c{Ta^Pi) - h 



M- 



Given any e > we increase Ri > Rq such that the second term on the right 
hand side is smaller than e. Next we choose R2 > 2Ri such that the first 
term is small. Now that Ri and R2 are fixed, the third term converges to 
zero by Lemma [3.61 and since Ta^p^ is minimizing the remainder follows suit. 
Therefore for i sufficiently large, 

/ Ta,p^dx = M-Mi,s>M-{-Chl,6o)-^e. 

The strong convergence of the potential energies now follows by Lemma 13.61 
By weak convergence po ^ a.e., and for any e>0 there exists R>0 such 
that 



M> / podx>M-e, 
Jbr 

in particular pQ L^CR."^) with J pQ = M. The functional p^ J'^{p)dx is 
convex, so by Mazur's lemma 2.13] and Fatou's lemma 



/ ^(po)dx<limsup / '^{Ta^pi)dx. 

J i—>oo J 
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Hence po G J\j with 



Wc(/)o) < limsupWc(Pi) = ^M) 

and the proof is complete. □ 

Remark, (a) Thm. 13.71 provides a minimizer pQ of the reduced energy- 
Casimir functional T-C^ under the assumptions (^1)~(^'3). By Thm. [T^ this 
minimizer can be lifted to a minimizer /q of the original energy-Casimir 
functional Tic ■ By Lemma 12.21 the function ^ satisfies the necessary as- 
sumptions if ^ which appears in the original Casimir functional satisfies the 
following ones: ^> G CH[0,oo[), $(0) = ^>'(0) = and 

(<I>1) <I> is strictly convex, 

($2) $(/) > Cf+^Z'' for / > large, 

($3) $(/) < C7/i+i/^'' for / > small, 

with growth rates k,k' £]0,1[. 

(b) As will be seen in the next section the mere fact that /o minimizes TCc is 
not sufficient for stability. However, let (/j) C J-'m be a minimizing sequence 
of Tic- By Thm. 12.31 (a) the sequence of induced spatial densities Pi = Pf^ 
is minimizing for Ti.Q. Choose a subsequence of {pi) (and (/j)) and shift 
vectors such that the assertions of Thm. 13.71 hold, and denote the shifted 
subsequence again by (/j). We claim that this sequence converges weakly to 
/q. Clearly, (/j) is bounded in L-^+^/'^(R^) with bounded kinetic energy, and 
Epot{fi) = EpotiPi)^Epot{po)- Any subsequence of (Ji) must therefore have 
a weakly convergent subsequence with weak limit /o which is a minimizer 
of Tic and induces the same spatial density po and potential Uq. But then 
by Thm. El /o = /o so that indeed fi^fo weakly in L1+^/'=(M^). 

(c) For A; > 1 one can still obtain stability results, cf. for the Kuzmin 
disk which corresponds to ^(/) = /^^^, i.e., k = k' = 2. However, the reduc- 
tion approach cannot work, because as we shall see in the last section this 
approach implies stability for the Euler-Poisson system where stability is 
probably lost at n = 2, i.e., k = l. 

4 Stability of minimizers 

Now that the existence of a minimizer is proven, we can explore its dynamical 
stability properties. So let po be as obtained in Thm. rTTI and /o as induced 
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by Thm. 12.31 A simple expansion shows that 

nc{f)-nc{fo) = d{f,fo) + E^ot{f-fo), (4.2) 
where for f£j-M and with the Lagrange multipher £^0 from Thm. [T^ (h). 

difjo) ■■= JJmf)-Hfo) + E{f-fo)]dvdx 

= JJ[^{f)-^fo) + {E-Eo){f-fo)]dvdx 

> ll[^'{fo) + {E-Eo)]{f-fo)dvdx>0 

with d{f,fo) = iff f = fo- For the positivity of d we use the strict con- 
vexity of ^ and the form of /o according to Thm. 12.31 (b); by that theo- 
rem the term in brackets vanishes on the support of /q. We recah that 
--E^pot(/) = (/0/,P/)pot = ||/0/||pot defines a norm on pf, cf. Lemma|TIl note 
that the right hand side in Eqn. (|4.2I) is (i(/,/o) — 1 1/9/ — /OqI Ipot- We obtain 
the following stability result; C^(M^) denotes the space of compactly sup- 
ported functions on M^. 

Theorem 4.1 Let /q be a minimizer of Tic on Tm obtained from a min- 
imizer po ofTiQ, and assume that the minimizer is unique. Then for any 
e>0 there exists 5>0 such that for any classical solution [0,T[b th^ f{t) of 
the fiat Vlasov-Poisson system with /(O) G C^(]R^) H^m and 1 1/(0)| li+i/^ = 
||/o||i+i/fc the estimate 

d{f{0)Jo) + \\Pfio)-Po\\lot<S 
implies that for each [0,T[ there exists a shift vector aGlR^ such that 

\\fit)-Tafo\\l+l/k + d{f{t),Tafo) + \\pf^t)-TaPo\\lot<^- 

Proof. Assume that the assertion were false. Then there exists e > 0, tj > 0, 
fj (0) G C2 n with 1 1 (0) 1 1 = 1 1 /o 1 1 such that for every j G N, 

d{fj (0) , /o) - i^pot if J (0) - /o) < J (4.3) 
but for any shift vector a G M^, 

II/,- (tj) - Tafo\\i+i/k + d{fj{tj),Tafo) - Epotifj (tj) - TJo) > e. (4.4) 
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Since Tic is preserved along solutions we have from ()4.2|) and ()4.3|1 that 

nc{Mtj))=nc{fjm^nc{fo), 

i.e., {fj{tj)) is a minimizing sequence. By Thm. and the remark at the 
end of the previous section there is a sequence of shift vectors (aj) C M? such 
that up to a subsequence, 

By (|4.2|1 this implies that d{Ta^fj{tj),fo)^0. For the convergence of 
||-||i+i/fc we use the fact that ||/j(t)||i+i/fc = ||/j(0)||i+i/fe = ||/o||i+i/fc for 
any t>0. By the remark, Tajfj{tj)^fo weakly in and hence 

TaJjitj)^fo strongly in L^+^/''{R^). But these convergence results for 
Tajfjltj) contradict dOJ. □ 

Remark, (a) The uniqueness assumption on /o in the above theorem is 
made mostly in order to avoid technical complications. It suffices if /o is 
isolated with respect to the topology of our stability estimate. If there should 
be a continuum of minimizers then the set of minimizers itself is stable; 
we refer to IHl Thm. 4] for such a formulation of the result in the three 
dimensional case. We are not aware of a case where there is a continuum 
of minimizers with fixed mass M. For a closely related variational problem 
it has been shown that the above stability estimate remains valid even then 

m- 

(b) As opposed to the three dimensional case |1U1I141[T3] there is no global ex- 
istence and uniqueness result to the initial value problem for the flat Vlasov- 
Poisson system yet. Hence our stability result is conditional in the sense that 
it holds as long as a suitable solution exists. A local existence and unique- 
ness result for smooth solutions with initial data in C^(M^) as well as a 
global existence result for weak solutions to the flat system was established 
in We could also carry out our stability analysis in the framework of 
these global weak solutions, but this would only bury the main ideas under 
technicalities. 

(c) By interpolation between and L^+^/^ we obtain a stability estimate 
for \\f{t) — Tafo\\p with p 1 -|- l//c]. If we assume that the initial pertur- 
bations have supports of uniformly bounded measure we can include the case 
p = l, if we assume a uniform bound on the L°° norm of the initial pertur- 
bations we can by interpolation include all [1 + 1/A;,cxd[; notice that both 
the measure of the support and the L°° norm are invariant under classical 
solutions of the Vlasov-Poisson system. 
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(d) The need for the shifts in the stabiUty estimate arises from the Gahlei 
invariance of the Vlasov-Poisson system. If /o is a steady state then for any 
fixed VgR'^ the function fo{x-tV,v-V) is a time dependent solution; /o 
is simply put into a uniformly moving coordinate system. But while the 
distance of this perturbation to the steady state grows linearly in t, it is 
arbitrarily close to the steady state at t = for V small. 



5 Connection to the Euler-Poisson system 

A self-gravitating matter distribution can be described on the microscopic, 
kinetic level represented by the Vlasov-Poisson system or on the macro- 
scopic, fluid level represented by the Euler-Poisson system. The reduction 
technique connects the stability problems for these two viewpoints. In the 
three dimensional situation this connection was observed in • In the flat 
case the corresponding Euler-Poisson system reads 

dp 

— + div{pu) = 0, 

du 

p— + (u-Vx)u = -VxP - pVxU, 
J \x-y\ 

with the equation of state 

Here p denotes the pressure of the fluid and u denotes its velocity field; 
the meaning of p and U is as before. If pQ is a minimizer of the reduced 
energy-Casimir functional TCq, then using the Euler-Lagrange identity in 
Thm. 12.31 (b) it is easy to check that po and the zero velocity field uo = 
solve the flat Euler-Poisson system. Clearly, the state {pq,uq) minimizes the 
energy 

n{p,u) = ^ j \u\^pdx + j ^'(p)dx + ^pot(p), 

among the states with J p = M. Formally, the energy is conserved along 
solutions of the Euler-Poisson system. An expansion about {po,uq) gives 

W(p,n)-W(/Jo,^io) = ^ j \u\'^P'^x + d{p,pi)) + E^ot{p- Pa), 
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where 

d(/>,po):= j [^ip)-^ipo) + iUo-Eo){p-po)]dx>0. 

Now the stabihty proof proceeds in the same way as in the Vlasov case. We 
can for every e > find a 6>0 such that for every solution 1 1— > {p{t),u{t)) of 
the flat Euler-Poisson system with p{0)(zJ^j, which preserves energy and 
mass, the initial estimate 

^ I \um'p(.0)dx + d{p{0),po) + MO)-po\\lot<S 

implies that as long as the solution exists and up to shifts in space, 

PO 1 1 l+l/n + 1 j W{t)\^p{t)dx + d{p{t),po) + \\p{t)-po\\l^,<e. 

Neither in the flat case nor in the three dimensional one is there an existence 
theory for global solutions of the Euler-Poisson system, which preserve all 
the necessary quantities, so the result is conditional in this sense. 
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